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EXPONENTIAL ERGODICITY OF THE JUMP-DIFFUSION CIR 

PROCESS 

PENG JIN, BARBARA RUDIGER, AND CHIRAZ TRABELSI 


Abstract. In this paper we study the jump-diffusion CIR process (shorted 
as JCIR), which is an extension of the classical CIR model. The jumps of the 
JCIR are introduced with the help of a pure-jump Levy process (Jt,t > 0). 
Under some suitable conditions on the Levy measure of ( Jt , t > 0), we derive 
a lower bound for the transition densities of the JCIR process. We also find 
some sufficient condition guaranteeing the existence of a Forster-Lyapunov 
function for the JCIR process, which allows us to prove its exponential er- 
godicity. 


1. Introduction 


The Cox-Ingersoll-Ross model (or CIR model) was introduced in [1] by John C. 
Cox, Jonathan E. Ingersoll and Stephen A. Ross in order to describe the random 
evolution of interest rates. The CIR model captures many features of the real 
world interest rates. In particular, the interest rate in the CIR model is non¬ 
negative and mean-reverting. Because of its vast applications in mathematical 
finance, some extenssions of the CIR model have been introduced and studied, see 
e.g. [3, 4, 13]. 

In this paper we study an extension of the CIR model including jumps, the 
so-called jump-diffusion CIR process (shorted as JCIR). The JCIR process is de¬ 
fined as the unique strong solution X := (X t ,t > 0) to the following stochastic 
differential equation 

dX t = a{6 - X t )dt + ay/X t dW t + dJ t , X 0 > 0, (1.1) 


where a, u > 0, 9 > 0 are constants, ( Wt,t > 0) is a 1-dimensional Brownian 
motion and ( J t ,t > 0) is an independent pure-jump Levy process with its Levy 
measure v concentrating on (0,oo) and satisfying 



(£ A \)u{d£ > ) < oo. 


( 1 . 2 ) 


The initial value Xq is assumed to be independent of (W t . t > 0) and ( J t ,t > 0). 
We assume that all the above processes are defined on some filtered probability 
space (fi, J 7 , (X) t >o,P). We remark that the existence and uniqueness of strong 
solutions to (1.1) is guaranteed by the main results of [6]. 
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The term a(9 — Xt) in (1.1) defines a mean reverting drift pulling the process 
towards its long-term value 6 with a speed of adjustment equal to a. Since the 
diffusion coefficient in the SDE (1.1) is degenerate at 0 and only positive jumps 
are allowed, the JCIR process ( X t ,t > 0) stays non-negative if X 0 > 0. This fact 
can be shown rigorously with the help of comparison theorems for SDEs, for more 
details we refer the readers to [6]. 

The JCIR defined in (1.1) includes the basic affine jump-diffusion (or BAJD) as 
a special case, in which the Levy process ( J tl t > 0) takes the form of a compound 
poisson process with exponentially distributed jumps. The BAJD was introduced 
by Duffie and Garleanu [3] to describe the dynamics of default intensity. It was 
also used in [4] and [12] as a short-rate model. Motivated by some applications in 
finance, the long-time behavior of the BAJD has been well studied. According to 
the main results of [8, 12], the BAJD possesses a unique invariant probability mea¬ 
sure, whose distributional properties was later investigated in [7, 9]. We remark 
that the results in [8, 9] are very general and hold for a large class of affine process 
with state space R + . The existence and some approximations of the transition 
densities of the BAJD can be found in [5]. A closed formula of the transition 
densities of the BAJD was recently derived in [7]. 

In this paper we are interested in two problems concerning the JCIR defined in 
(1.1). The first one is the transition density estimates of the JCIR. Our first main 
result of this paper is the derivation of a lower bound on the transition densities 
of the JCIR. Our idea to establish the lower bound of the transition densities is 
as follows. Like the BAJD, the JCIR is also an affine processes in R + . Roughly 
speaking, affine processes are Markov processes for which the logarithm of the 
characteristic function of the process is affine with respect to the initial state. 
Affine processes on the canonical state space R™ x R n have been investigated 
in [2, 4, 10, 11]. Based on the exponential-affine structure of the JCIR, we are 
able to compute its characteristic function explicitly. This enables us further to 
represent the distribution of the JCIR as the convolution of two distributions. The 
first distribution is known and coincides with the distribution of the CIR model. 
However, the second distribution is more complicated. We will give some sufficient 
condition such that the second distribution is singular at the point 0. In this way 
we derive a lower bound estimate of the transition densities of the JCIR. 

The other problem we consider here is the exponential ergodicity of the JCIR. 
According to the main results of [8] (see also [12]), the JCIR has a unique invariant 
probability measure 7r, given that some integrability condition on the levy measure 
of (Jt,t > 0) is satisfied. Under some sharper assumptions we show in this paper 
that the convergence of the law of the JCIR process to its invariant probability 
measure under the total variation norm is exponentially fast, which is called the 
exponential ergodicity. Our method is the same as in [7], namely we show the 
existence of a Forster-Lyapunov function and then apply the general framework of 
[14, 15, 16] to get the exponential ergodicity. 

The remainder of this paper is organized as follows. In Section 2 we collect 
some key facts on the JCIR and in particular derive its characteristic function. In 
Section 3 we study the characteristic function of the JCIR and prove a lower bound 



EXPONENTIAL ERGODICITY OF THE JUMP-DIFFUSION CIR PROCESS 


3 


of its transition densities. In Section 4 we show the existence of a Forster-Lyapunov 
function and the exponential ergodicity for the JCIR. 


2. Preliminaries 


In this section we use the exponential-affine structure of the JCIR process to 
derive its characteristic functions. 

We recall that the JCIR process (X t ,t > 0) is defined to be the solution to 
(1.1) and it depends obviously on its initial value Xq. From now on we denote by 
{Xf,t > 0) the JCIR process started from a initial point x > 0, namely 

dX? =a{0-X?)dt + cjy/xfdW t +dJ t , X% = x. (2.1) 


Since the JCIR process is an affine process, the corresponding characteristic 
functions of {X£,t > 0) is of exponential-affine form: 

E[e uX * ] = u G U := {u G C : SRw < 0}, (2.2) 

where the functions <j>(t, u) and u) in turn are given as solutions of the gener¬ 
alized Riccati equations 


I d t <j>(t,u) = F(i/>(t,u)), <j>(0,u) = 0, 

| dtip(t,u) = R(ip(t,u)), ip{ 0> u ) = u 

with the functions F and R given by 

F(u) =adu+ f (e M? — 1 )u(d£), 

J (0,oo) 

2 2 

o/ \ ° u 

n[u) -—-- au. 


Solving the system (2.3) give <j)(t,u) and ip(t,u) in explicit form: 


(2.3) 


and 




ue at 

1 - f|«(l - e-»‘) 


(2.4) 


2 aO / cr 2 

1 ~ 2a 111 — e 


))+/ l (e^ {s ’ u) -l)v(dOds. (2.5) 

JO J (0,oo) V ' 


Here the complex-valued logarithmic function log(-) is understood to be its main 
branch defined on C — {0}. According to (2.2), (2.4) and (2.5) the characteristic 
functions of {X£,t > 0) is given by 


E[i 


■uXf 




• exp 


(t 


xue 


- f ^ 1 - e 


• exp j J ^ e ^( s <“) — l^u(d£)ds^ 


for any u £U. Here the complex-valued power function z := 
is also understood to be its main branch defined on C — {0}. 


—at 


) 


( 2 . 6 ) 


exp(-^rlogz) 
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We denote by 


/y^ _ 2 ad 

I :=(l - —u(l — e~ at )) • exp 


( _* 

\ 1 — — 

^ 2a 


sue 


t(l - e~ at ) 


II:=exp(jf (e^ s ’ u ) - l)i/(df)ds) (2.7) 

the first and the second line of the characteristic function (2.6) respectively. 

According to the parameters JCIR process we look at two special cases: 

2.1. Special case i): v = 0, no jumps. Notice that v = 0 just yields the clas¬ 
sical CIR model (Y t , t > 0) satisfying the following stochastic differential equation 

dY t x = a{6-Y x )dt + a^fdW t , Y x =x> 0. (2.8) 

It follows from (2.6) that the characteristic function of ( Y f x ,t > 0) coincides with 
I, namely for u £U 


E[e uY *} =(1 - ^-u( 1 - e~ at )) & • exp 
=1 


(r 


- h u ( l - e ~ at ) 


(2.9) 


It is well known that the classical CIR model ( Y x ,t > 0) has transition density 
functions f{t,x,y) given by 


f{t, x, y) = ne “ vI q (2(uv)i) 


( 2 . 10 ) 


for t > 0, x > 0 and y > 0, where 

2a 


cr 2 (l - e“ ot ) ’ 


v =ny, 


2 ad , 

Q = a 1) 


and / 9 (-) is the modified Bessel function of the first kind of order q. For x = 0 the 
formula of the density function p(t, x, y) is given by 




( 2 . 11 ) 


for t > 0 and y > 0. 


2.2. Special case ii): 0=0 and x = 0. We denote by (Z t ,t > 0) the JCIR 
process given by 

dZ t = —aZ t dt + a\/~ZtdWt + dJt, Zq = 0. 

In this particular case the characteristic functions of ( Z t ,t > 0) is equal to II, 
namely for u £ Id 
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where 


4>(t,u) = 


1 - ! k u ( 1 ~ e at ) 

One can notice that II resembles the characteristic function of a compound poisson 
distribution. 


3. A lower bound for the transition densities of JCIR 

In this section we will find some conditions on the Levy measure v of (Jt, t > 0) 
such that an explicit lower bound for the transition densities of the JCIR process 
given in (2.1) can be derived. As a first step we note that the law of Xf,t > 0, 
in (2.1) is absolutely continuous with respect to the Lebesgue measure and thus 
possesses a density function. 

Lemma 3.1. Consider the JCIR process (Xf,t > 0) (started from x>0) that is 
defined in (2.1). Then for any t > 0 and x > 0 the law of Xf is absolutely contin¬ 
uous with respect to the Lebesgue measure and thus possesses a density function 
P{t,x,y), y> 0. 

Proof. As shown in the previous section, it holds 

E[e uX *} = I 11= E[e uY * X ) ■ E[e uZt ], 

therefore the law of Xf. denoted by pxf , is the convolution of the laws of Y t x and 
Z t . Since ( Y f x ,t > 0) is the well-known CIR process and has transition density 
functions f(t, x, y), t > 0, x, y > 0 with respect to the Lebesgue measure, thus Hx? 
is also absolutely continuous with respect to the Lebesgue measure and possesses 
a density function. 

□ 

In order to get a lower bound for the transition densities of the JCIR process 
we need the following lemma. 

Lemma 3.2. Suppose that f^ Q ^ £ ln(i)u(d£) < oo. Then II defined in (2.7) 
is the characteristic function of a compound Poisson distribution. In particular, 
P(Z t = 0) > 0 for all t. > 0. 

Proof. Since 



If we rewrite 


( fe as u \ ( au \ 

eXP l 1 _ f d.(l_e-^)J =eXP l^J’ 


(3.1) 
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where 


2 at; 


T 2 (e as — 1) 


2ae a 


> 0 


> 0 , 


(3.2) 


' ' a 2 (e as — 1) 

then we recognize that the right-hand side of (3.1) is the characteristic function of 
a Bessel distribution with parameters a and /?. Recall that a probability measure 
Ma,/3 on (R + ,£>(R + )) is called a Bessel distribution with parameters a and /3 if 


Va,p(dx) = e a 5o(dx) + (3e a — • Ii(2y/a/3x)dx, (3.3) 

where So is the Dirac measure at the origin and I\ is the modified Bessel function 
of the first kind, namely 


oo ( \ 2\ * 

Ur) = -Y ) 

u k'-ik + iy.’ 


r G 


For more properties of Bessel distributions we refer the readers to section 3 of [7]. 
Let p, a< p denote the characteristic function of the Bessel distribution with 
parameters a and /3 which are defined in (3.2). It follows from [7, Lemma 3.1] 
that 

£e _os u 


Mo 


e („) = exp(- § 2IL)= e xp( I -^ 


(1 - e~ as J 


:)• 


Therefore 


A = 


n (/V/3 (u) - 1 )v{d£)ds 

v 0,oo) 
rt 


0 J (0,oo' 


( 


e p-u — e 


- 1 


^is(d(;)ds. 


Set 


A := I f (l — e a )is(dt;)ds 

J 0 J (0,oo) 
rt 


0 J (0,oo) 


1 — e o- 2 (e a,s -l) 


^u(d£)ds. 


(3.4) 


If A < oo, then 


A = 


g/J-U _ g 


/0 J (0,oo 

rt 


= A 


A 


0 J (0,oo) 


^ v(d£,)ds — A 
(g^=« — e~ a )v(dt;)ds — 1^. 


The fact that A < oo will be shown later in this proof. Next we show that the term 
— / / (e^ 1 — e~ a )is(dl;)ds can be viewed as the characteristic function of 

* Jo J(O.oo) 
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a probability measure p , which we now construct as mixtures of the following 
measures 

a 


m. 


t ,p{dx) := fde “ /3x w— • Ii(2^/a(3x)dx, x > 0, 


where I\ is the modified Bessel function of the first kind. Noticing that the measure 
m a< p is the absolute continuous component of the measure p, a< p in (3.3), we easily 
get 

m a ,p(u) = p. a ,p{u) - e~ a = - e ~ a , 

where rh a ,p{u) '■= / 0 °° e ux m a ^{dx) for u €U. Since a and /? are functions in £ 
and s, we can define a measure p on R + as follows: 

P := \ f f rn a ,p v(d£)ds. 

A do 0,oo) 

By the definition of the constant A in (3.4) we get 


1 

A 

1 

A 

= 1 


[ [ rn at/3 (R + )v(d£)ds 

J 0 J (0,oo) 

f j (1 — e~ a )is(dt;)ds 

J 0 J (0,00) 


i.e. p is a probability measure on R + and for u € U 

p{u) = f e ux p(dx) 

J (0,oo) 


1 


A 


(e^“ — e a )v(dt;)ds. 

/0 J (0,oo) 

thus A = A (p(u) — 1) and E[e uZt ] = j s characteristic function of a 

compound Poisson distribution. 

Now we verify that A < 00 . Since 

rt 


A = I f (l — e Q )u(d£)ds 

J 0 J (0,oo) 
rt 


'0 *d(0,oo) 
r t 


I( 


2 ag 


' (0,oo) J 0 


We introduce the change of variables 


l — e <r' j (e a *-t)'ji/(d£)ds 


2 ag 


l—e a^(e as — 1) 

2a£ 


t 2(e as — 1) 


:= y and then get 


dy = -- 


2a£ 


cr 2 (e as — l) 2 


2 cr % 2 < , n , 

= -y b7l— + !)ds. 
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Therefore 


2 a£ 


A = 


f (0,oo) J oo 

/ vm r (i-e-y) 

J (0,oo) J „, 2 ?5 


(l-O o-... , — dy 


2a^y + a 2 y 2 

25 


2( e at_i) 
OO 


'(o, 

Cr 2 ( e at _ ]\ 

where 6 := ---. Define 


[ vmTil-e-y) 

J CO, ool J 4 


2 a£y + cr 2 y 2 

25 


dy 


2 a£y + cry 


■dy, 


2a 


/»oo 

M(C):= y £ (l-e-*); 


2 ? 


rdy, 


' 2a5y + <J 2 y 2 

then M( 5 ) is continuous on ( 0 , oo) and as 5 —>• 0 we get 
Af(0 = / e (l-e"") 1 ^ 


2a5y + (j 2 i / 2 


/ OO 

(1 - e-*) 


2a5y + a 2 y 2 


< 2 ^ 


< 25 


y 


§ 2a£y + a 2 y 2 


dy + 25 


1 


1 


'« 2a5 + o- 2 y 


dy + 25 


1 2 a 5 y + a 2 y 2 
0 1 


dy 


1 2 a 5 j/ + <r 2 y : 


rdy. 


For y £ [§, 1] we have 


25 


4 2a5 + cr 2 y 


1 7 2 ? 

- dy = ^ 


G L 

25 


ln(2a5 + cr 2 y) 

^ 2 5 


n 1 


ln(2a5 H-?-) — ln(2a5 + o' 2 ) 

Cl- 0 

= ci 5 + c 2 5 ln(^) < c 3 5 ln(i) 


and as 5 —t 00 , 


/*00 

M(5) < y s (l-e-«) 


25 


25 


fi 2 aty + G 2 y 2 


2a£y + a 2 y' 
dy < 25 


-,dy 


rvh-^i 

■ 1 


- 2 /- 


We yield 


< 


25 <5 2d 

= ■ 7 = -9 := c 4 < 00 . 

5 cr 


A < c 3 / 5ln(i)u(d5) +c 4 [ 

Jo S Jl 


i <j 2 y 2 


lu(d5) < 00. 
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□ 


With the help of the previous lemma we can easily prove the following propo¬ 
sition. 


Proposition 3.3. Let p(t,x,y), t > 0,x, y > 0 denote the transition densities of 
the JCIR process (Xf,t > 0) defined in (2.1). Suppose that J ^ £ln(|)u(d£) < 
oo, then for all t > 0, x, y > 0 we have 

P(t,x,y) > C(t)f(t,x,y), 

where C(t) > 0 for all t > 0 and f(t,x,y) are transition densities of the CIR 
process (without jumps). 


Proof. According to Lennna (3.2), we have P(Zt = 0) > 0. Define 

C(t) :=P(Z t = 0). 

Since 

E[e uX *} =111 = E[e uY *} ■ E[e uZt ], 

the law of Xf , denoted by px*, is the convolution of the laws of Y t x and Z t . Thus 
for all A G B(R+) 


Pxf (A) = 


(A - y)p Zt ( dy) 


> 


'{ 0 } 


PY? (A - y)pz t (dy) 


> pyx (A - {0 })pz t ({0}) 

> C(t) ■ py*(A) 

> C(t) [ f(t,x,y)dy , 

J A 

where f(t, x, y) are the transition densities of the classical CIR process given in 
(2.10). Since A G B(K+) is arbitrary, we yield 


P(t,x,y) > C(t)f(t, x, y) 


for all t > 0, x, y > 0. 


□ 


4. Exponential ergodicity of JCIR 

In this section we find some sufficient conditions such that the JCIR process 
is exponentially ergodic. Our method relies on the lower bound of the transition 
densities we have derived in the previous section and the existence of a Forster- 
Lyapunov function. 

Lemma 4.1. Suppose that ^ £u(g?£) < oo. Then the function V(x) = x, x > 0, 
is a Forster-Lyapunov function for the JCIR process defined in (2.1), in the sense 
that for all t > 0, x > 0, 

E[V(Xf)\ < e~ at V(x) + M, 
where 0 < M < oo is a constant. 
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Proof. Since jixf = /ly* * Hz t , thus 

E{X?] = E[Y*] + E{Z t ]. 

Since (Y t x ,t > 0) is the CIR process starting from x, it is known that /iy.® is a 
non-central Chi-squared distribution and thus E\Y f x ] < oo. Next we show that 
E[Z t \ < oo. 

Let u £ (—oo,0). Then by using Fatou’s Lemma we get 


E[Z t \ = £[lim 


e uZt - 1, 


u—>0 U 

uZ t _ i E\e uZt 1 — 1 

< lim inf M--] = lim inf ^-- 

u —>0 U u —^0 U 


Remember that 


= e A(u) , 


-<( 

du\ 


a t o _ y ue _ 

J (e 1 -^ (1 - e_as) “ - 1 )j/(d^)ds) 

£r«“(l - £ • (1 - e"“» - Zue~ as (£ • (1 - e~ as )) 


then we have for all u < 0 

£ue~ as 


"fsr ( 1_e as ' )u —1)) = 


)) 


(! - h ' (1 -e~ as )u) 2 


< 


?e" 


i 1 ~ f-(f - e ~ as )uy 


<£e~ 


and further 


/ 0 J (0,oo) 

Thus A (u) is differentiable in u and 

rt 


f f £e as v(d^)ds < oo. 

J 0 J (0,oo) 


A'(0) = [ [ te~ as v(dOds = -—— [ 
Jo 0,oo) a J(( 


(0,oo) 




It follows that 


E[Z t \ < lim inf 


II(u) - 11(0) 


= II'( 0) = e A(0) • A'(0) 

1 - e~ at f 

= -/ MdO- 

^ J (0,oo) 


Therefore under the assumption / £z/(c££) < oo we have proved that E[Z t \ < 

J (0,oo) 


oo. Furthermore 


E[z ,i = |.( £|e -. ] ) 


1 — e 


u=0 


' (0,c 


&{<%)■ 
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On the other hand 


E[e^]=( 1 —^(l-e-))--.exp(- ^ 


with a similar argument as above we get 

E[Y t x ] = A = 0(1 - e" at ) + xe" at . 

Altogether we get 

E[X X ] = E[Y t x ] + E[Z t ] 

1 — p~ a ^ 

= (1 — e~ at )(6 +-) + xe~ at 

a 

< 9+-+xe~ at , 

a 


namely 


E[V(X?)] <e + -+ e~ at V(x). 


□ 


Remark 4.2. If ^ £i/(d£) < oo, then it follows from [ 8 , Theorem 3.16] and [12, 
Proposition 3.1] the existence and uniqueness of an invariant probability measure 
7 i for the JCIR process. 

Let || • ||tv denote the total-variation norm for signed measures on R + , namely 

IImIItu = sup {|m(A)|}. 

AeB(R+) 

Let P t (x,-) := P(X'f G •) be the distribution of the JCIR process at time t 
started from the initial point x > 0. Now we prove the main result of this paper. 

Theorem 4.3. Assume that 

[ £ v {d£) < 00 an d f £ln(—)i/(d£) < oo. 

4(1,oo) 4(0,1) ? 

Let 7 r be the unique invariant probability measure for the JCIR process. Then the 
JCIR process is exponentially ergodic, namely there exist constants 0 < (3 < 1 and 
0 < B < oo such that 

\\P t {x, •) — tt\\tv < B(x + l)/3 4 , t > 0, xgR+. (4.1) 

Proof. Basically, we follow the proof of [16, Theorem 6.1]. For any 8 > 0 we 
consider the ^-skeleton chain rff := Xf 5 , n € Z+. Then (?)“)„ e z + is a Markov 
chain on the state space R+ with transition kernel P s (x,-) and starting point 
?)q = x. It is easy to see that the measure tt is also an invariant probability 
measure for the chain (?^) ra€ z + , x > 0 . 

Let V(x) = x, x > 0. It follows from the Markov property and Lemma 4.1 that 
E[VK + iM, V x i,---,V x n }= [ V(y)P 5 (r 1 x ,dy)<e- aS V( V x ) + M, 

J M_l_ 
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where M is a positive constant. If we set Vo := V and V n := n £ N, then 

E[Vi] <e~ aS V 0 (x) + M 

and 

E[Vn+i\Vo,Vi, •••,<]< e~ a5 V n + M, n £ N. 

Now we proceed to show that the chain {rj^)nei + i x > 0, is A-irreducible, strong 
aperiodic, and all compact subsets of R + are petite for the chain (^)„ e z + . 

“A-irreducibility”: We show that the Lebesgue measure A on R + is an irre- 
ducibility measure for (r)%)ne z+- Let A £ ,B(R+) and A (A) > 0, then it follows 
from Proposition 3.3 that 

P[Vi e A\r% =x\= P(Xf € A) > C7(5) [ f(S,x,y)dy > 0, 

J A 

since /(<5, x,y) >0 for any x £ R + and y > 0. This shows that the chain (?^) n ez+ 
is irreducible with A being an irreducibility measure. 

“Strong aperiodicity” (see [14, p.561] for a definition): To show the strong 
aperiodicity of (?^) raS z 0 , w e need to find a set B £ B(M+), a probability measure 
m with m(B) = 1, and e > 0 such that 

L(x,B)> 0, x £ 1 + (4.2) 

and 

P(rji £ A) > e - m(A), x £ C, A £ B(R + ), (4.3) 

where L{x,B) := P{rj * £ B for some n £ N). To this end set B := [0,1] and 
g(y) := inf^.g^ i ]f(6,x,y), y > 0. Since for fixed y > 0 the function f(5,x,y) 
strictly positive and continuous in x £ [ 0 , 1 ], thus we have g(y) > 0 and 0 < 
f(o,i] 9{y)dy < 1. Define 

m(A) := - - 1 f g(y)dy , A £ S(R+). 

J( 0 , 1 ] 9(y)dy J A n( 0 , 1 ] 

Then for any x £ [0,1] and A £ B(R + ) we get 
P(vi e A) =p(Xf £ A) 

>C(6) [ f(5,x,y)dy>C(8) [ g(y)dy = C(6)m(A) f g{y)dy , 
JA JAn(0,l] J(0,1] 

so (4.3) holds with e := C(8) f^ 0 1 , g(y)dy. 

Obviously 

L(x, [0,1]) > P(if, £ [0,1]) = P(X% £ [0,1]) > C(8) [ f(S , s, y)dy > 0 

d[0,l] 

for all x £ R + , which verifies (4.2). 

“Compact subsets are petite”: We have shown that A is an irreducibility mea¬ 
sure for {rin)nez + - According to [14, Theorem 3.4(h)], to show that all compact 
sets are petit, it suffices to prove the Feller property of (? 7 ^)ngz + , x > 0, but this 
follows from the fact that (? 7 ^)nez + is a skeleton chain of the JCIR process, which 
is an affine process and possess the Feller property. 
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According to [14, Theorem 6.3] (see also the proof of [14, Theorem 6.1]), 
the probability measure 7r is the only invariant probability measure of the chain 
(Vn)ne z+, x > 0, and there exist constants /3 G (0,1) and C G (0, oo) such that 

\\P Sn (x, •) — tt\\tv < C(x + l)(3 n , iiG Z+, iGl+ 

Then for the rest of the proof we can proceed as in [16, p.536] and get the 
inequality (4.1). □ 
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